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We present Integral Field Unit GMOS-IFU data of the compact HII galaxy 
UM 408, obtained at Gemini South telescope, in order to derive the spatial dis- 
tribution of emission lines and line ratios, kinematics, plasma parameters, and 
oxygen abundances as well the integrated properties over an area of 3"x4".4 
equivalent with ~750 x 1100 pc located in the central part of the galaxy. The 
starburst in this area is resolved into two giant regions of about 1".5 and 1" (~375 
and ~250 pc) diameter, respectively and separated 1.5-2" (~500 pc). The ex- 
tinction distribution concentrate its highest values close but not coincident with 
the maxima of Ha emission around each one of the detected regions. This indi- 
cates that the dust has been displaced from the exciting clusters by the action of 
their stellar winds. The ages of these two regions, estimated using H/? equivalent 
widths, suggest that they are coeval events of ~5 Myr with stellar masses of 
~10'^Mq. We have also used [0111]/H/3 and [Sll]/Ha ratio maps to explore the 
excitation mechanisms in this galaxy. Comparing the data points with theoreti- 
cal diagnostic models, we found that all of them are consistent with excitation by 
photoionization by massive stars. The Ha emission line was used to measure the 
radial velocity and velocity dispersion. The heliocentric radial velocity shows an 
apparent systemic motion where the east part of the galaxy is blueshifted, while 
the west part is redshifted, with a relative motion of ~ 10 km s^^. The veloc- 
ity dispersion map shows supersonic values typical for extragalactic Hll regions. 
Oxygen abundances were calculated from the [Olll]AA4959,5007/[OlIl]A4363 ra- 
tios. We derived an integrated oxygen abundance of 12+log(0/H)=7.87 summing 
over all spaxels in our field of view. An average value of 12-|-log(0/H)=7.77 
and a difference of A(O/H)=0.47 between the minimum and maximum values 
(7.58±0.06-8.05±0.04) were found, considering all data points where the oxygen 
abundance was measured. The spatial distribution of oxygen abundance does 
not show any significant gradient across the galaxy. On the other hand, the bulk 
of data points are lying in a region of ±2a dispersion (with a=0.1 dex) around 
the average value, confirming that this compact Hll galaxy as other previously 
studied dwarf irregular galaxies is chemically homogeneous. Therefore, the new 
metals processed and injected by the current star formation episode are possibly 
not observed and reside in the hot gas phase, whereas the metals from previous 
events are well mixed and homogeneously distributed through the whole extent 
of the galaxy. 

Subject headings: galaxies: individual (UM 408) — galaxies: abundances — 
galaxies: dwarf — galaxies: ISM. 
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Introduction 



HII galaxies (or Blue Compact Dwarfs, BCDs) are low-mass and metal-poor galaxies 
(I/5O-I/3Z0), experiencing strong episodes of star forr aation, characterized by the pres- 



ence of bright emission lin es on a faint blue con t inuum ( Terlevich et al 



991 



Kehrig et al. 



Cairos et al. 



2003 



2004f). Several studi es fe.g. lPapaderos et al.lll996l : lTelles fc Terlevich]|l997 
Westera et al.l 120041 ) have shown the existence of an old stellar population underlying the 
present burst in most of these galaxies, indicating that these objects are not young systems 
forming their first generation of stars. 

The star formation activity in HII galaxies is concentrated in several luminous star 
clusters (~l-30 pc in size) spread over the galaxies. Some of these clusters have been asso- 
ciated with supe r star clusters (SSCs), similar to those found in int eracting galaxies, such 



as the antennae (IWhitmore fc Schweizei 



galax i es such as Henize 2-1 (IConti fc Vacca 



1995 



I99J and NGC 1569 (e.g., lO'Connell et al. 



Whitmore et al.lll999l ). and some local dwarf 



I994J : lArp fc Sandagd Il985l ). During their evolution, these clusters ionize the interstellar 
medium (ISM) causing the formation of giant HII regions (GHIIRs), and release a consider- 
able quantity of me chanical energy into the ISM, characterized by supersonic motions (e.g, 
Melnick et al.lll987l ). sweeping up the surrounding medium, removing the gas and dust from 
the star formation site, while producing structures such as bubbles and supershells. Ulti- 
mately, the evolution of these young clusters also causes the freshly produced metals to be 
dumped into the ISM. 

Oxygen and other a elements are produced by massive stars (>8-10 Mq) and are re- 
leased into the ISM during their supernova phase. These metals will be dispersed in the 
whole galaxy and m ixed via hydrodynamic mechanisms in time scales of few 10*^ yrs (e.g., 
Tenorio-Tagld Il996l ) . Since the spatial distribution of abundances of O, N, etc depends on 
the recycling time of the ISM, the spatial variation of these abundances give important 
insights about these processes. The spatial distribution of emission line ratios and of the 
abundance of certain element s in dwarf galaxies hav e been studied by different authors (e.g., 
Kobulnicky fc Skillmaru Il996l . 119971 : iLee et al.l |2006| ) . In these works the radial and spatial 
distribution of the 0/H abundance, obtained from the analysis of their H II regions, have 
been used as a chemical evolution indicator. Localized nitrogen self enrichment has been 
measured in a few cases and attributed to the action of strong w inds produced by Wolf-Rayet 
(W-R) stars fe.g.. IWalsh fc Rovlll989l : iKobulnickv et al.lll997h. _ b ut no oxygen localized en- 
richment systems have been confirmed (IKobulnicky fc Skillmaru Il996l ). In GHIIRs, where 
only one star formation region is present, it is usually assumed that the abundance of the 
galaxy and their ISM is well represented by the abundance of this single region, assuming 
that the metals are well mixed and distributed across the galaxy. 
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The details of the structure of HII galaxies are important to understand their ionization 
mechanisms, star formation feedback and chemical enrichment. These issues have been ac- 
tively addressed in s tudies of the brightest and nearest systems . The large and he terogeneous 
HII galaxy sample (ILoose fc ThuanI Il986l : iKunth et al.l Il988l : iTelles et al.l 119971 ) includes a 
significant number of far less studied objects at larger distances, with small apparent sizes 
and with relatively low surface bri ghtness. Typically , HII galaxies contain one or a few star 
formation regions, but H/3 images fiLagos et al.ll2007n of some compact objects have revealed 
that the central region probably hosts a myriad of unresolved star clusters. The present fa- 
cilities, with high spatial resolution and high efficiency instruments on large telescopes (8m), 
allow us to study these more distant objects to derive their observed properties and relate 
them to the better known nearby galaxies, thus having a handle on the intrinsic properties 
of star formation as well as possible evolutionary effects. 

In this paper we use integral field spectroscopy with Gemini Multi-object Spectrograph 
and the Integral Field Unit (GMOS-IFU) at Gemini South in order to study the spatial 
distribution of emission lines, their ratios, extinction, abundances and the kinematic prop- 
erties of the gas in the interstellar medium of the compact dwarf galaxy UM 408. UM 408 
belongs to a subset sample of ra re HII galaxies with low metal abu ndances (i.e. Z< 1/20 
Zq) with 12 + \og{0/H) = 7.66 (IMasegosa et al.lll994j : lTelleslll995l ). compact morphology, 
and very small effe ctive radius (Rsf/) of on ly 2".l from previ ous morphology and surface 
brightness studies (ITelles fc TerlevichI 119971 ). On other hand, IPustilnik et al.l ( 120021 ) using 
high res olution spectroscop y derived an abundance of 12 + log[0/H) = 7.93 and HI obser- 
vat ions dSalzer et al.ll2002l ) reveal an HI mass of log(M///)=8.815 Mp. UM 408 is cataloged 
bv ISalzeir(ll989h as a Dwarf HII Hot Spot with a redshift of 0.0120 (v=3598 km s"^) with 
coordinates RA= 02^^ 11"^ 23%, DEC= 02° 20' 30"(J2000) at a distance of ~ 46 Mpc, and V 
apparent magnitude of 17.38. With the present observations, we note that the central part 
of UM 408 encompasses two main star formation regions (here named A and B), as shown 
by the g-band acquisition images in Figure [H Another faint region, called C, can be seen in 
the outer parts of the galaxy. These regions were all unresolved in previous studies. 

The observations and data reduction procedures are discussed in Section [2l Section [3] 
gives the results obtained from the 2D emission line maps. In Section H] we discuss our results 
and our conclusions are presented in section O 



^Value obtained from NED 
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UM 408 

z = 0.012 




Fig. 1. — g-band acquisition image of the galaxy UM 408. The image is 11". 5x11". 5. 
The rectangle indicates the field of view of 3"x4".4 (~750 x 1100 pc) used in this work 
to construct the monochromatic maps. The two most prominent regions (A and B) are 
indicated in the image. The faint region (C in the plot) located ~3".0 to the East from B is 
placed outside the field of view (1":^^ 250 pc). 



2. Observations and Data Reductions 



The observat ions were performe d with the Gemini Multi-Ob ject Spectrograph CMOS 
faook et al.ll2004j ) and the IFU unit jAUington-Smith et al.lbooi hereafter GMOS-IFU) at 
Gemini South Telescope in Chile during the night of August 21 and December 31 2004, using 
the grating B600+_G5323 (B600) and on January 01, 2005 with the grating R600+_G5324 
(R600) in one slit mode, covering a total spectral range from ~3021 to 7225 A. The GMOS- 
IFU in one slit mode composes a pattern of 750 hexagonal elements, each with a projected 
diameter of 0".2, covering a total 3". 5 x5" field of view, where 250 of these elements are 
dedicated to sky observation. The detector is formed by three 2048x4608 CCDs with 13.5 
Hm pixels, with a scale of 0.073"/pixel. The CCDs create a mosaic of 6144x4608 pixels with 
a small gap between the chips of 37 columns. 



Table 1 shows the observing log, indicating the instrumental setup, airmass and exposure 
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Fig. 2. — ID blue and red parts of the spectrum, obtained as the sum of all fibers, in the 
rest-frame wavelength. In this Figure, only the lines used in the analysis are indicated. 



times, the dispersion, the final resolution and the seeing (FWHM) of each observation. 
The data were reduced using the Gemini package version 1.8 inside IRAfi. All science 
exposures, comparison lamps and spectroscopic twilight and GCAL flats were overscan/bias 
subtracted and trimmed. The spectroscopic GCAL fiats were processed by removing the 
calibration unit plus GMOS spectral response and the calibration unit uneven illumination. 
Twilight fiats were used to correct for the illumination pattern in the GCAL lamp flat 
using the task gfresponse in the GMOS package. The twilight spectra were divided by 
the response map obtained from the lamp fiats and the resulting spectra were averaged in 
the dispersion direction, giving the ratio of sky to lamp response for each fiber. The final 
response maps were obtained then by multiplying the GCAL lamp fiat by the derived ratio. 
The resulting extracted spectra were then wavelength calibrated, corrected by the relative 
fiber throughputs, and extracted. The residual values in the wavelength solution for 60-70 
points using a 4*^ or S*'^ Chebyshev polynomial typically yielded rms values of ~0.10 A 
and ~0.15 A for both gratings respectively. The final spectra cover a wavelength interval 
of ~302 1-5823 A and ~4371-7225 A for data taken with the B600 and R600 gratings, 
respectively. The final spectral resolutions and dispersions are shown in Table 1 (columns 5 



^IRAF is distributed by NO AO, which is operated by the Association of Universities for Research in 
Astronomy Inc., under cooperative agreement with the National Science Foundation. 
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Fig. 3. — Subsets from the spectra in Figure 2 showing details of selected emission lines 
considered in this study. 

and 6). 

The flux calibration was performed using the sensitive function derived from observation 
of the star LTT 3864 for both gratings in December 31 and January 1. Standard stars were 
not observed in August. Therefore, the blue part of the spectra was constructed using only 
the observation obtained in December 2004. The 2D data images were transformed into 3D 
data cube (x,y,A) using the gfcube routine and resampled as square pixels with 0".l of spatial 
resolution. The emission line fluxes were measured using the IRAF task fitprofs by fitting 
Gaussian profiles. From these, we then created the emission line maps used in our analysis 
(see Section [3]). 
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At shorter wavelength the IFU spectroscopy and long-sht observations suffer a spatial 
translation produced by differential atmospheric refraction (DAR). This effect is wavelength 
dependent, and is produced by the deviation of the light when it passes through the at- 
mosphere due to the variation of the air density as a function of elevation . In order to 



obtain fluxes corrected of PAR, we used a similar proce dure as described in lArribas et al. 



( I1999I ) and also applied by IWestmoquette et al.l (l2007al ). First we splitted the data cube 
in monochromatic images, and for each image we calculated the centroid of an unresolved 
point source in the field of view, that in our case corresponds to region A. Finally, the various 
monochromatic maps were aligned by shifting the centroids to the same position. 

In Figure [2] we show the resultant spectrum for the two gratings used, obtained from the 
sum of all GMOS-IFU fibers. Figure [3] shows a detailed view of the selected emission lines 
considered in this study. The spectral resolution of the GMOS-IFU observations allow us to 
resolve the [OII]AA3726,29 doublet lines in a limited number of apertures with a AA=2.32 A 
between the peak of the lines in the integrated spectrum (see Figure [3]). The low intensity of 
[OII]AA3726,29 lines and [OIII]A4363 are the largest sources of uncertainties in our results. 
An automatic procedure to measure these emission line fluxes was not possible, therefore 
these lines were measured manually using splot. 



3. Results 

3.1. Emission line, continuum and EW(H/9) maps 

We used a flux measurement procedure described in the previous section to produce the 
maps for the following emission lines: [SII]A6731, [SII]A6717, [NII]A6584, Ha, [OIII]A5007, 
[OIII]A4959, H/5, [OIII]A4363, H7, [OII]AA3726,29, Ha and R/3 continuum and EW(H/3), 
in a rectangular field of 3"x4".4 (see Figured]). In order to reduce the number of data 
points without loosing spatial information and S/N, the data were resampled to 0".2 in 
spatial resolution. Typical S/N ratios for individual pixels (spaxels in our data cube) for 
[OIII]A5007 are >700 and > 500 for regions A and B, respectively. The [OIII]A4363 line 
shows a typical S/N ratio >8 and >4 for regions A and B, respectively. The signal to noise 
for each emission line is given by the ratio between the emission line peak and the RMS of the 
adjacent continuum. Figure H] shows some of the observed emission line maps. All emission 
and continuum maps in this Figure are in units of x 10~^* erg cm~^ s~^, except for the H/3 
equivalent width map, which is in units of A. Here 0".2 ^ 50 pc (distances in this work were 
computed assuming a Hubble constant of Ho=72 km s~^ Mpc~^). Since most of the emission 
lines were measured using an automatic procedure, we filtered the maps assigning the value 
erg cm~^s~^ to all pixels with S/N<3. 
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Fig. 4.— Observed maps for [SII]A6717, [NII]A6584, Ha, [OIII]A5007, [OIII]A4363, H7 and 
[OII]AA3726,29 emission lines, H/3 continuum and EW(H/3). Fluxes in units of xlO~^^ erg 
cm -2 s"^ and EW(H/3) in units of A. Ha emission line contours are overplotted on all maps. 
The maximum Ha emission is placed over region A and is indicated in the maps by an X 
symbol. The regions in the [OIII]A5007 map mark the areas considered as regions A and B. 
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The maps in Figure H] show essentially the same spatial distribution, presenting two well 
defined and regular regions, that we have labeled A and B. The maximum peak intensity 
of the Ha emission line in region A is indicated in the maps by an X symbol. The recom- 
bination lines Ha, H/5 and H7 have very similar spatial distribution to the forbidden lines 
in Figure HI but there is a slight difference in the spatial distribution of [OIII]A4363 and 
[OII]AA3726,29 over region B. The last two images in Figure H] corresponds to the spatial 
distribution of the H/3 continuum and the EW(H/3), respectively. It is clear from the maps 
that the continuum peak is not coincident with the maximum in Ha emission. This has been 



previously reported by o ther authors (e.g.. lSargent &: Filippenkdll99ll : iMaiz-Apellaniz et al. 



19981 : iLagos et al.l 120071) in some galaxies with single and multiple s t ar for ming regions. In 



NGC 4214, ISargent fc Filippenkol fll99lh and iMaiz-Apellaniz et al.l fll998f ) found that the 
continuum peak emission was displaced in relation to the line emitting regions where a 
strong Wolf-Rayet (WR) feature (Hell A4686 A) was also seen. No WR signatures were 
detected in the spectra of UM 408. 

The largest EW(H/?) values in the galaxy are associated with region A, indicating that 
this region is younger than region B. On the other hand, the B.f5 continuum map (see Figure 
S]) shows the largest values over region B. In Section HJ the EW(H/3) values are used as 
a proxy for the starburst ages in order to derive the integrated physical properties of the 
galaxy and their star forming regions. T he error associate d to ea ch emission line (cTj) was 
calculated using the relationship given in ICastellanos et al.l (120021 ) : 



EW 

7VA' 



where is the standard deviation in the local continuum associated to each emission line, 
is the number of pixels, EW is the equivalent width of the line and A is the instrumental 
dispersion in A (see Table 1 ). These estimated errors will be used later in the derivation of 
uncertainties in the oxygen abundance determination and integrated properties. 

Finally, given that the seeing was 0.8-0.9" (~4 pixels in the emission line maps), each 
aperture or pixel of 0.2" should be correlated with its neighbor, and the pixel-to-pixel varia- 
tions in the maps produced from the observed emission line are likely due to the uncertainties 
attached to the measurement and do not represent real variations. Therefore, we smoothed 
the emission line ratio, temperature and oxygen abundance maps using a 3x3 box. 
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Fig. 5. — Extinction coefficient distribution, obtained from the emission lines Ha/H/3 ratios. 
The regions enclosed by circles mark the location of other extinction peaks, probably related 
to other unresolved star clusters. The maximum Ha emission is placed over region A and is 
indicated by an X symbol. Ha contours are overplotted on the maps. 

3.2. Extinction 

The logarithmic reddening parameter c(H/?) was calculated from the ratio Ha/H/3, 
where the intrinsic value of 2.87 for case B recombination was assumed. Thus, the cor- 
rected emission lines are calculated as 



/(A) F(A) 



liHP) F{H/3) 



(2) 



where 1(A) and F(A) are the dereddened flux and obs erved fl u x at a given wavelength, 
respectively and f(A) is the reddening function given by ISeatoru (119791 ) as parametrized by 
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Howarth I (Il983l ). using appropriated values for the Milky Way (IKobulnicky et al.l 119971 ) 



Figure O shows the smoothed spatial distribution of c(H/3) for the emission line ratios 
mentioned above. The reddening parameters c(H/3) were set to 0.0 for unrealistic values of 
Ha/H/? < 2.87. We observe a range of extinction of ~0.02 to 1.04 with average c(H/3) value 
of 0.76 and a standard deviation of 0.17. The lowest extinction values are located in the outer 
part of the field of view (see Figure [5]). We note also that the maximum extinction c(H/5) is 
clo se to, but not coincident wit h the peak of Ha e r nission . A similar result has been reported 
by iMaiz-Apellaniz et al.l ( 19981 ) and ICairos et al.l ( 20021 ) for the star formation knots in the 
galaxies NGC 4214 and Mrk 370, respectively. Apparently, the current starburst episode in 
this galaxy is sweeping the gas and dust out of the center into the surrounding regions. The 
emission line distribution as well as the c(H/3) distribution over region A and its elongated 
shape suggest that this region may contain two or more star forming complexes instead of 
only one. These complexes are indicated with circles in the extinction maps. 

Finally, as the extinction across the galaxy does not show a uniform distribution, con- 
centrating the highest values around the star forming regions, the fluxes of each pixel in the 
emission line maps were corrected using their corresponding c(H/5) value. 



3.3. Nebular emission line diagnostic diagrams 



The standard BPT diagnostic diagrams (IBaldwin et al.lll98ll ) have been used to ana- 



lyze the possible excitation mechanisms present in UM 408. Figure [6] (a, b and c) shows 
the following emission line ratio maps: [OllI]A5007/H/3 ([01I1]/H/?), [SII]AA6717,6731/Ha 
([SII]/Ha) and [NII]A6584/Ha ([Nil] /Ha). Figure E] shows that [0III]/H/3 and [Sll]/Ha 
ratios have essentially the same spatial distribution with inverse trends. The largest values 
of [0111] /H/9 and the lowest values of [SII]/Ha are seen in the east part of the galaxy, near 
region B. The emission line ratios measured within these regions are common to high exci- 
tation Hll regions with values ranging from log([OIII]/H/3)=0.45 to 0.80 (all map), 0.65 to 
0.76 (region A) and 0.72 to 0.80 (region B), and from log([SII]/Ha)=-1.21 to -0.62 (all map), 
-1.17 to -0.85 (region A) and -1.21 to -1.01 (region B). Additionally, the low ionization ratios 
[Nil] /Ha and [SII]/Ha increase from the center outwards. 

Figure [7] shows the log[OIII]/H/3 vs. log[SII]/Ha diagnostic diagram, identifying three 
different regions in UM 408 using the original data points (non-smoothed). Red and green 
crosses correspond to regions A and B, respectively, and the black crosses correspond to the 
outer parts of the galaxy. The integrated value that represents the sum over all spaxels in 
our field of view is given by the blue star symbol. The solid line in Figure [7] divides the 
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Fig. 6.— Emission line ratios: (a) Log [OIII]A5007/H/?, (b) Log [SII]AA6717,6731/Ha and 
(c) Log [NII]A6584/Ha;. The maximum Ha emission is placed over region A and is indicated 
in the maps by an X symbol. Ha contours are overplotted on Figure a). 



diagram into two regions that correspond to line ratios explained solely by photoionization 
by massive stars and to line ratios where an additi onal source of excitatio r i mus t be present, 
as in the case of Active Galactic Nuclei (AGNs; lOsterbrock fc Ferlandl l2006l ). From the 
distribution of the data in the diagnostic diagram we conclude that the HII regions in UM 
408 are produced by photoionization by massive stars. 



3.4. Density, temperature and oxygen abundance 



The first step in the abundance derivation is to obtain the electron density and temper- 
ature. The electron density is obtained using the emission line ratio [SII]A6717/[SII]A6731, 
and the electron temperature from the ratio [OIII]AA4959,5007/[OIII]A4363. In Figure [S^ 
we can see that th e ratio [SH]AA6717,6731 is typically greater than 1 which indicates a 
low density regime (lOsterbrock fc Ferlandl l2006l ) . hence, we assumed an electron density of 
Ug ~100 cm~^ for all apertures in our calculations. To obtain the electro n temperature 
Te(OIII) we used the five level atomic model FIVEL (jPe Robertis et al.lll987l ). implemented 
under the IRAF STS package NEBULAR. This program has been developed to obtain the 
physical conditions in a low-density nebula, given appropriate emission line ratios. Figure 
IHb shows the [OIII]AA4959,5007/[OIII]A4363 emission line ratio used to calculate the elec- 
tron temperature Te(OIII), and Figure [HI:; shows the resulting spatial distribution of electron 



-14- 




Log [SII]XX6717,6731/Ha 



Fig. 7. — [OIII]/H/? versus [SII]/Ha diagnostic diagram for the original data points presented 
in Figure Ell) and c). Red crosses correspond to region A, green crosses to region B and 
the black crosses to the outer part of the galaxy. The integrated value for the galaxy is 
represe nted by the blue star. The solid line is the dividing line between HII galaxies and 
AGNs JOsterbrock fc Ferlandl 120061 ). 



temperature Te(OIII). The range of valid data points varies from 1.20 to 1.70 x 10^ K with 
a standard deviation of 0.1 x 10^ K for the original data set and from ~1.37 to ~1.53 x 10^ 
K with a standard deviation of 0.03 x 10^ K for the smoothed data points. The total oxygen 
abundances for each aperture are obtained assuming the contribution of O"^ and O"*"^ ions, 
thus we have that 



O 
H 



0+ o- 



H+ 



(3) 



where the Tp( OII) temperature a nd the O"^ and 
sions given by IPagel et al.l (1l992l ) . 



ions are obtained assuming the expres- 



Figure [9|i shows the spatial distribution of the oxygen abundances over the regions 
where the emission line [OIII]A4363 was detected. The oxygen abundance values in units of 
12+log(0/H) range from 7.70±0.08 to 7.84±0.04 with a maximum error of ~0.12 dex, an 
average value of 7.77 and a standard deviation of 0.10 dex. However, the smoothed map 
has values from ~7.69 to ~7.86 with a standard deviation of ~0.04 dex. Figure [9b shows 



- 15 - 



□ ) SIIX6717/SII)>6751 b) [OIII]X4959,5007/[OIII]X4363 c) Te[01II] x 10" K 




arcsec arcsec arcsec 



Fig. 8. — (a) [SII]A6717/[SII]A6731 ratio map used to calculate the electron density, (b) 
[OIII]AA4959,5007/[OIII]A4363 emission line ratio, which is used to calculate the electron 
temperature Te(OIII). (c) Te(OIII) electron temperature in units of lO^K. The maximum 
Ha emission is placed over region A and is indicated in the maps by an X symbol. Ha 
contours are overplotted on the maps. 

the radial distribution of oxygen abundance with respect to the Ha peak emission within a 
region < 0.4 kpc. Statistically the bulk of the original data points are lying in a region of 
2(T dispersion (ct=0.1 dex) around the average value. The 2ct dispersion is represented by 
the dotted lines. A least fit square with a slope of -0.14±0.06 dex kpc~^ was found, using 
all original data points in Figure [9]d. The error in the gradient is obtained directly from the 
linear regression. In the same Figure we show the radial distribution of the smoothed data 
points as light blue triangles. These results indicate that there is no significant variations 
across the galaxy. At most, we can say that there is a very marginal gradient of a decreasing 
abundance from the center outwards, indicating that, on average, the highest abundance 
values are found near the peak Ha emission. 



3.5. Velocity field 

The internal structure of HII galaxies and GHIIRs, as viewed in the ionized gas, is 
directly associated with mechanisms of photoionization, magnetic field induced turbulence, 
and feedback by the current episode of massive star formation and evolution, i.e. radiative 
shocks, stellar and supernovae driven-winds. All of it, under the influence of the gravi- 
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b) Radial distribution 
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Fig. 9. — Oxygen abundances [12+Log(0/H)]: (a) Spatial distribution. The isocontours 
display the Ha emission. The maximum Ha emission is placed over region A and is indi- 
cated in the map by an X symbol. Ha contours are overplotted on the maps, (b) Radial 
distribution. The central dotted line represents the 12+log(0/H) average value of 7.77, and 
the solid line the integrated value of 12-|-log(0/H)=7.81 in the region where the emission 
lines [OIII]A4363 were measured. The difference between the maximum and minimum values 
in the original data set (black dots) is equal to 0.47 dex. Light blue triangles represents the 
smoothed 12-|-Log(0/H) radial distribution. A least square fit using the original data points 
is indicated by the red line. 



tational potential of the complex of gas and stars. The presence of expanding structures 
(shells, loops and bubbles) is very common in GH IIRs (e.g., the HST image on NGC 604, 



Tenorio-Tagle et al.l 12000 : 



Chu fc KennicuttI Il994l . for 30 Doradus). These structures have 



also been observed in the prototypical HIT galaxies as II Zw 40 and their resulting nested 
filaments have bee n explained by the effects of photoionization and stellar winds in an in- 
homogeneous ISM (ITenorio-Tagle et al.ll2006l ). UM 408 shows, in the Ha emission line map 
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[FigureH]) and in the acquisition broad-band image (Figur e [T]), an outer regular shape with 



no signs of large-scale disruption or a galactic wind (see IWestmoquette et al.l l2007bl . and 
references therein for the case of NGC 1569). The monochromatic maps, presented in this 
study, do not provide clear evidence of expanding structures, shell like features or filaments. 

In order to study the impact of the star formation on its internal kinematic properties, 
we derived the spatial distribution of radial velocity v,, (heliocentric) and velocity dispersion 
cr (FWHM/2.355) by fitting a single Gaussian to the Ho; emission line profiles. Figure [TD] 
shows the radial velocity derived from the shifts of the Ha line peak, and dispersion velocity 
maps calculated by using the following relation 



a\Ha) = crlf,^ - a^^^t - a^^, (4) 

where aobs is the observed Ha velocity dispersion, ainst is the instrumental dispersion, and ath 
is a correction for thermal broadening [^JkTe/mH ~9.1 km s~^ for hydrogen at a temperature 
of Te=10^ K). We have used a value of 33.4±1.0 km s~^ for the instrumental dispersion 
by comparing the nominal value (from Table 13. ip with that measured directly from the 
calibration lamps in the grim R600 near Ha, and assigning the error as the rms of the 
average. 

The velocity field (Figure [TOh) shows an apparent systemic motion where the east part 
of the galaxy (near knot B) is blueshifted, while the west part (near knot A) is redshifted, 
with a relative motion of ~ 10 km s~^. This motion may contribute to broaden the integrated 
line profile over the extent of the whole galaxy. However, it is not enough to explain the 
total integrated supersonic line width observed of a = 18.6 ± 1.5 km s~^. Figure [TOb shows 
the velocity dispersion (line width) map over the observed field, ranging values from ~ 10 
km s~^ to ~ 30 km s~^. The velocity dispersion varies little across the field, but some points 
are worth noting. This map is characterized by a subsonic area near region A, while the 
highest values are located in the outer part of the field of view near region B. The one pixel 
line width over knot A (the brightest region) has a low a of 17.1±1.0 km s~^. This value is 
identical to a synthetic aperture of 2.7" over knot A. This syn t hetic aperture was chosen to 



mimic the single fiber high dispersion observation of iBordald (120041 ) who find (T=19.6±0.4 
km s~^, that is consistent with our results. 
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a) V,(Ha) km s'' 



b) a(Ha) km s'' 
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Fig. 10. — (a) Ha radial velocity in units of km s~^. (b) Ha Velocity dispersion a. The 
value showed in this figure is obtained considering an instrumental amst and thermal ath 
correction over the observed value aobs- The maximum Ha emission is placed over region A 
and is indicated in the maps by an X symbol. 



4. Discussion 

4.1. The integrated properties of UM 408 

We have measured the fiux in two different apertures enclosing more than ~60% of the 
observed Ha flux of the galaxy in order to study the properties of the star formation regions 
detected in this galaxy. One aperture for region A (48% of the measured Ha flux) and another 
for region B (13% of the measured Ha flux). These regions have a diameters of about 1".5 
and 1" that correspond to approximately 375 pc and 250 pc, respectively. Regions A and B 
are indicated in the [OIH]A5007 map of FigureHl These regions are formed by the data points 
where the [Oin]A4363 emission line is detected. To o btain the integrated properties of the 



galax y and compare with previous observations (e.g. iMasegosa et al.lll994 



Pustilnik et al. 



20021 ) we derived the integrated spectra summing over all pixels in our data cube. Table [2] 
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summarizes the resultant emission line values, obtained using the different synthetic aper- 
tures, including the total [OIII]A4363 aperture over the galaxy (in regions where [OIII]A4363 
emission line is detected and measured) . These meas urements are in reaso r iable agreement 



with those derived with lo ng slit spectroscopy (e.g., iTerlevich et al.lll99ll : iMasegosa et al. 
1994 : iPustilnik et al.ll2002l ) considering aperture differences and observational uncertainties. 
A larger di screpancy however occ urs between our observed [OII]AA3726,29/H/3 value of 0.52 
and that of IPustilnik et al.l (120021 ) of 1.45. However, their value must be overestimated since, 
from their Figure 1, we clearly see that H/? > [OII]A3727. 

The observed EW(H/3) values in UM 408 are ~102 A and ~57 A for regions A and B 
and ~67 A for the integrated galaxy, respectively. The values for regions A and B have been 
calculated by taking the average over all pixels in each region. Using these EW(H/9) values we 
can esti mate the cluste r ages by comparison with the values obtained from STARBURST99 
model (ILeitherer et al.l 119991 ). The resulting ages assuming an instantaneous burst and a 
metallicity of Z=0.004 for a Salpeter IMF with a mass limit of IM© to lOOM© are 4.83 Myrs 
and 5.00 Myrs for regions A and B, respectively. These ages are consistent with both regions 
being formed simultaneously. We found integrated luminosities (assuming a distance of 
D=46.8 Mpc; Pustilnik et al. 2002) of Log L(Ha)=39.15, 38.58 and 39.47 erg s"^ for regions 
A and B and for the integrated galaxy, respectively. These correspond to a star formation 
rate of SFR=0.011, 0.003 and 0.024 Mq yr~^, and with a number of ionizing photons log 
Q(H'^)=51.01, 50.45 and 51.34 photons s~^ for regions A , B and for the integrated galaxy, 
obtained from the relationships given by iKennicuttI (|l998l ). The masses of these regions can 
be estimated by scaling the number of ionizing photons with the models predicted for the 
corresponding ages. The stellar masses are 5.17 x 10^ M© and 2.29 x 10^ M0 for regions 
A and B, respectively. Table [3] shows the derived integrated properties of the galaxy and its 
star forming regions. 

Since its discovery, UM 408 was classified as a low abun dance galaxy with an i ntegrated 
value of 12+Log(0/H)=7.63 for low resolution spectro scopy (jMasegosa et al. 1994) and 7.93 
for high resolution spectroscopy (IPustilnik et al.ll2002l ). In this work, we found an integrated 
abundance (summing over all spaxels in our field of view) of 12+Log(O/H)=7.87±0.05. 
The synthetic apertures previously mentioned, yield an integrated oxygen abundance of 
12-|-Log(O/H)=7.82±0.05 and 7.84±0.06 for regions A and B, respectively. Measuring the 
abundance in an aperture equivalent with the area where the emission line [OIII]A4363 was 
measured we obtain a value of 7.81±0.05, the same value obtained in the peak of the Ha 
emission (see Figure [Hb). This result reflects the fact that the integrated values are light 
weighted and dominated by the galaxy peak emission. 



Finally, despite the fact that a detailed kinematic analysis is beyond the scope of the 
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present paper, our results confirm that the core dominates the kinematic inforn aation of the 
star forming region, and present a low a value, as found by other works (e.g., iTelles et al. 



2OOII . and references therein). Therefore, a single Gaussian fit to the line profile of any 



aperture encompassing the brightest knot will measure a representative line width of the 
dominant supersonic motions, somewhat deconvolved by the effects of stellar and SN me- 
chanical energy input. The effects of stellar or SN feedback will contribute to the broadening 
of the line profiles, but will only be detectable in the lower density regions at low intensi- 
ties. One means to recognize these effects is through the i dentification of inclin e d ban ds in 
the diagnos t ic dia grams of intensity vs. a as proposed by iMunoz-Tunon et al.l (119961 ) and 



Yang et al.l (119961 ). However, the compactness of UM 408 and the small range of reliable 



line widths make the use of the diagnostic diagrams difficult to interpret in this case. The 
Ha line profile is symmetric and well represented by a single Gaussian, and does not show 
prominent low intensity wings in either apertures. An attempt to deblend a possible addi- 
tional component to this line profile produced only an upper limit of abroad < 100 km s~^ 
with F{a)broadlF{a) 

narrow 

<3%. 



4.2. Distribution of the oxygen abundance and their comparison with other 

dwarf galaxies 

The uniform behavior of 0/H abundance in scales of hundreds of pc in low mass galaxies 
as in UM 408 is not without precedent, and is comparable with the observed variation in 
some dwarf galaxies by other authors. In NGC 4214 (Dwarf Irregular/ Wolf- Rayet galaxy) 
Kobulnicky fc Skillmaru (119961 ) showed that there is no localized O, N or He abundance 
gradi ents, but they found, nea r the y oungest region, an oxygen overabundance of 0.095±0.019 
dex. iKobulnicky fc SkillmanI (Il997l ). using long slit spectroscopy to study the interstellar 
medium of the dwarf irregular galaxy NG C 1569, fa i led to find evidence for 0/H gradient 
from the recent star formation activity. 



Lee et al.l (120061 ) studied the local group dwarf 



irregular galaxy NGC 6822 and measured a difference of A(O/H)=0.53 dex between the 
maximum (8.43±0.2) and minimum value (7.90±0.1), with a mean oxygen abundance of 
12+log(O/H)=8.11±0.1 using only the HII regions where the [OIII]A4363 emission lines 
were detected. A gradient of -0.14±0.07 dex kpc~^ was measured by Lee et al. in a radius 
<1.4 kpc, and a slope of -0.16 dex kpc^^ was obtained if measurements from the literature 
were included, sho wing the existence of a possible radial gradient in oxygen abundance. 
Izotov et al] J2006[ l using VLT/GIRAFFE in the ARGUS mode, observed the dwarf galaxy 
SBS 0335-052E, one of the most metal po or galaxies with an integrated oxygen abundance of 
12-|-log(O/H)=7.30 (IMelnick et al.lll992l ). The spatial distribution of the oxygen abundance 
in the ISM of SBS 0335-052E vary from 7.00±0.08 to 7.42±0.02, representing a difference 
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of A(O/H)=0.42 dex. The maximum value of oxygen abundances in this galaxy does not 
correspond to the position of one of the identified SSCs, but Izotov et al. find a slight trend 
of a decreasing abundance, interpreted as possible self-enrichment. However, they argue that 
the errors in the calculated oxygen abundances and Te, Ug, etc, not considering instrumental 
and observational uncertainties, lead to consider the oxygen abundance variations in SBS 
0335-052E rn a y not be statistically significant. More recently, using IFU-PMAS observations. 



Kehrig et al.l (120081 ) studied the spatial distribution of some metals over the galaxy IIZw70. 
In this study they find a difference of A(O/H)=0.40 dex between the maximum (8.05±0.06) 
and minimum (7.65±0.06) values of oxygen abundance. In the case of UM 408, the bulk of 
our observed data points are distributed in a region of ±2cr dispersion around the average 
value with a difference between the maximum and minimum values of 0.47 dex. This is, 
somewhat, in good agreement with previous results in other dwarf galaxies, using in some 
cases different techniques. We note that the gradient found in our work, calculated in a 
spatial scale of hundreds of pc, is similar to the gradients calculated in other studies in scales 
of kpc. Finally, the two giant regions A and B show a difference of oxygen abundance of 
only A(O/H)=0.02 dex, which indicates that these regions show identical chemical properties 
within the errors. Given the ages (~5 Myr) and stellar masses (~1O*^M0) of these regions, 
we expect that hundreds of SNs have exploded, producing eventually localized enrichments. 
But the absence of chemical overabundances in the ISM of UM 408 and in the dwarf galaxies 
studied in the literature lead to conclude that the population of young clusters have not 
produce preferentially localized overabundances. 

These results are compatible with the interpretation that the newly synthesized metals 
from the current star formation episode may not be in the warm phase of the ISM, and thus 
are not observed at optical wavelengths. The metals that are observed, however, ought to 
be from previous events, and may be well mixed and homogeneously distributed over the 
whole extent of this low mass and compact galaxy. Note that in UM 408, there is no bar 
induced rotat ion, or shear, to pro duce the metal dispersal and mixing as in more massive 



galaxies (e.g.. iRoy &: Kuntlj|l995l ). and thus there must be another mechanism responsible 
for the large-scale dispersal and mixing with the ISM. One such hydrodynamic mechanism 
for the transport and mixing of the metals produced by c ompact bursts of star formation 



into the ISM has been proposed by iTenorio-Tagld (119961 ). This model predicts that the 



energy injected by core-collapse supernovae change the physical conditions (density and 
temperature) of the interstellar medium, while undergoing a long excursion into the galactic 
volume. The injected matter is first thermalized, near the starburst, as it goes through a 
reverse shock. This generates the giant pressure that allows for the built up of superbubbles, 
kpc scale structures that may grow for up to 50 Myr. Afterwards, once the SN II phase is 
over, the giant superbubble gas begins to cool down by radiation. This happens first within 
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the parcels of gas that present the largest densities, promoting their thermal instability. 
The sudden loss of pressure within densest parcels of gas would immediately trigger the 
appearance of re-pressurizing shocks emanating from the low density hot gas. The process 
leads to a plethora of dense condensations that inevitably will fall and settle within the main 
body of the galaxy. The process of cooling and dispersal onto the main body of the galaxy 
will occur within a few 10® yr. For a complete mixing, the model assumes that a further star 
formation episode has to occur and through photoioniz ation the metal enrich condensations 
will expand and be fully mixed with the ISM (see also IStasinska et al.l 120071 1. This scenario 
seems plausible and compatible with our analysis for the case of UM 408. 



5. Conclusions 

We used GMOS-IFU spectrum of the compact HII galaxy UM 408, in order to derive 
the spatial distribution of emission line ratios, extinction c(H/3), radial velocity v^, velocity 
dispersion a, electron temperature and oxygen abundance (0/H) as well the integrated 
properties over an extended region of 3"x4".4 equivalent to 750 x 1100 pc in the central 
part of the galaxy. Below, we summarize our results and conclusions. 

1. The observed region of the galaxy includes two giant HII regions not resolved in pre- 
vious studied, here labelled A and B. The sizes of these ionized regions are ~375 and 
~250 pc for regions A and B, respectively. Another region labelled C was found to the 
East of region B. We do not observe large scale structures as in other GHIIRs and HII 
galaxies in the monochromatic maps. 

2. The distribution of dust content in the galaxy was derived from the Balmer line decre- 
ment (Hq;/H/3). The c(H/5) distribution concentrates the highest values close but not 
coincident with the peak of Ha emission in each ionized cluster (A and B). The dust 
seems to be displaced from the star formation regions by the action of the star cluster 
winds. 

3. We used [0III]/H/9 and [SII]/Ha ratios to investigate the possible excitation mechanism 
over the ISM of the galaxy. Comp aring the data points i n the diagnostic diagram 



with the theoretical model given by lOsterbrock &: Ferlandl (120061 ) . we found that all 



measured Ha flux arises from gas photoionized by massive stars 

The velocity field (v^) shows an apparent systemic motion where the east part of the 
galaxy (near knot B) is blueshifted, while the west part (near knot A) is redshifted, with 
a relative motion of ~ 10 km s~^ and a difference between the maximum a minimum 
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values of ~33 km s ^. The velocity dispersion map shows supersonic values, typical 
for extragalactic HII regions, ranging values from ~ 10 km s^^ to ~ 30 km s~^. 

5. The ages of the two giant regions detected in this study, were estimated from their 
EW(H/5) and the STARBURST99 models, suggesting that they are coeval events of 
~5 Myr with stellar masses of ~10'^ Mq. We see a marginal difference in the integrated 
oxygen abundance between these regions of A(O/H)=0.02 dex. 

6. Finally, as found in other nearby dwarf galaxies, we do not observe a gradient in oxygen 
abundance across the compact galaxy UM 408. The bulk of the observed data points 
are lying in a region of ±2cr dispersion {a=0.1 dex); therefore, the new metals formed 
in the current star formation episodes are not observed and reside probably in the hot 
gas phase (T~10^ K), whereas the metals from previous star formation events are well 
mixed and homogeneously distributed through the whole extent of the galaxy. 

All results presented here are suggestive that UM 408 is an unevolved low metallicity 
dwarf galaxy, undergoing a simultaneous episode of star formation over the whole optically 
observed extension, and it is a genuine example of the simplest starbursts occurring in 
galactic scale, possibly mimicking the properties one expects for young galaxies at high 
redshift. 
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Table 1. Observational setup 



Grating 


Central Wavelength 


Airmass 


Exposure Time 


Resolution 


Dispersion 


Seeing 




(A) 




(seconds) 


(A) 


(A/pixel) 


(") 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


B600 


4300 


1.24 


1 X 1800 


1.98 


0.45 


0.8 


R600 


5850 


1.47 


2x2000 


1.67 


0.47 


0.9 
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Table 2. Observed emission line fluxes for the integrated galaxy and regions A and B . 

Fluxes normalized to F(H/9). 





Integrated galaxy 


[OIII]A4363 measured 


Region A 


Region B 




(1) 


(2) 


(3) 


(4) 


[SI I] A6731 


0.24±0.01 


0.20±0.03 


0.20±0.03 


0.20±0.04 


[SII] A6717 


0.32±0.01 


0.28±0.04 


0.27±0.03 


0.27±0.04 


[Nil] A6584 


0.13±0.01 


0.11±0.03 


0.11±0.02 


0.10±0.03 


Ha 


5.94±0.07 


6.07±0.29 


6.14±0.28 


5.87±0.28 


[Oil I] A5007 


5.24±0.21 


5.33±0.30 


5.30±0.25 


5.60±0.36 


[OIII] A4959 


1.71±0.18 


1.72±0.11 


1.71±0.11 


1.80±0.11 


[OIII] A4363 


0.06±0.01 


0.07±0.01 


0.07±0.01 


0.07±0.02 


H7 


0.29±0.01 


0.31±0.04 


0.31±0.04 


0.36±0.05 


[Oil] AA3726,29 


0.52±0.17 


0.30±0.05 


0.32±0.04 


0.27±0.05 


F(H/3) 


39.70±0.20 


24.93±0.86 


18.63 ±0.62 


5.09±0.17 


Note. — Fluxes 


normalized to F(H/3) 


in units of xlO~^® erj 


— 1 —9 

5 s cm . 


Column (1): 



Integrated galaxy correspond with the values obtained from the sum of all spaxels in the data 
cube. Column (2): [OIII]A4363 measured correspond with sum of all apertures in the maps 
where [OIIIJA4363 was measured. Columns (3) and (4): Regions defined in Figure |4] 
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Table 3. Integrated properties of the observed area of UM 408. 





Integrated galaxy 


[OIII]A4363 measured 


Region A 


Region B 




(1) 


(2) 


(3) 


(4) 


c(H/3) 


0.93±0.02 


0.96±0.06 


0.97±0.06 


0.91±0.06 


Te{OIII) X lO* (K) 


1.40±0.08 


1.45±0.14 


1.45±0.12 


1.44±0.16 


Te(OII) X lO"* (K) 


1.32±0.11 


1.34±0.18 


1.34±0.16 


1.34±0.21 


He (cm~3) 


~100 


~100 


~100 


~100 


log[OIII] A5007/H/3 


0.69±0.02 


0.70 ±0.02 


0.69±0.02 


0.72±0.03 


[SI I] /Ha 


0.09±0.01 


0.08±0.01 


0.07±0.01 


0.08±0.02 


O+/H+xl05 


1.16±0.26 


0.64±0.09 


0.70±0.09 


0.57±0.10 


O++/H+X105 


6.35±0.58 


5.88±0.68 


5.85±0.60 


6.31±0.82 


12+log(0/H) 


7.87±0.05 


7.81±0.05 


7.82±0.05 


7.84±0.06 


log L(Ha) (erg s"-"-) 


39.47 


39.27 


39.15 


38.58 


log (number of photons) (photons s^^)^ 


51.34 


51.14 


51.01 


50.45 


SFR(Mq yr~i)t 


0.024 


0.015 


0.011 


0.003 


EW(H/3) (A) 


62 


89 


102 


57 


Age (Myr)tt 




4.73 


4.83 


5.00 



Note. — t Number of photons and SFR were obtained using the relations given by iKennicut^ lll998l '). tt Ages 
were obtained using the STARBURST99 predictions model, assuming an instantaneous star formation and a Salpeter 
IMF from 1 to 100 Mq . 



